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FOREWORD 


This  report  was  prepared  by  the  Mapping  and  Charting  Research  Laboratory 
of  the  Ohio  State  University  Research  Foundation  under  USAF  Contract  No.  AF 
IS  (600)-90.  The  contract  is  administered  by  the  Mapping  and  Charting  Branch 
of  the  Photographic  Reconnaissance  Laboratory,  Weapons  Components  Division, 
Wright  Air  Development  Center,  Wrights  Patters on  Air  Force  Base,  Ohio.  Mr. 

A.S.  Rosing  is  Project  Engineer  on  the  project  applicable  to  the  subject  of 
this  report. 

Research  and  Development  Order  Nos.  R-6S5-I44,  "Research  in  Photo gramme try 
and  Geodesy  for  Aeronautical  Charting,"  and  R-6S3-5S,  "Aeronautical  Charting 
Systems,"  are  applicable  to  this  report. 

This  report  was  originally  initiated  at  the  Ohio  State  University  Research 
Foundation  as  OSURF  Technical  Paper  No.  161. 
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ABSTRACT 


The  form  of  a  single  level  surface  that  envelops  all  attracting  matter 
and  the  value  of  the  potential  or  of  the  gravitational  attraction  at  a  definite 
point  on  or  outside  the  surface  determine  uniquely  the  field  of  force  on  and 
outside  of  the  surface.  This  statement  holds  good  even  when  we  add  to  the 
gravitational  attraction  the  centrifugal  "force"  due  to  uniform  rotation.  Hence, 
if  we  assume  that  the  surface  of  the  earth  is  an  exact  ellipsoid  of  rotation 
under  its  own  attraction  combined  with  the  centrifugal  force  of  uniform  rotation 
about  the  minor  axis,  the  field  of  force  on  and  outside  this  ellipsoid  is  unique¬ 
ly  determined,  if  we  assume,  for  instance,  the  value  of  gravity  at  the  equator. 

This  determination  is  carried  out  by  means  of  a  rather  special  set  of  curvi¬ 
linear  coordinates.  In  terms  of  these  coordinates  the  value  of  the  potential 
at  any  point  in  exterior  space  and  the  value  of  gravity  are  expressible  in 
closed  form  in  terms  of  Lengendre  functions  of  the  second  kind  with  Imaginary 
argument,  and  these  again  are  expressible  in  terms  of  the  elementary  functions. 

In  practice,  however,  it  is  found  more  convenient  to  expand  the  Lengendre 
functions  in  series,  also  the  formulas  dependent  on  them.  See  formula  (32)p.l2 
for  the  potential  expressed  in  closed  form,  and  formulas  (45),p.l 5*  and 
formulas  (6l)  and  (63),p.20  also  formulas  (?6)p.23  for  the  values  of  gravity 
on  the  ellipsoid. 

Various  incidental  developments  are  given.  This  is  first  of  a  series  of  at 
least  three  papers  on  the  general  subject. 


PUBLICATION  REVIEW 

The  publication  of  this  report  does  not  constitute  approval  by  the  Air 
Force  of  the  findings  or  conclusions  contained  therein.  It  is  published  only 
for  the  exchange  and  stimulation  of  ideas. 
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THE  GRAVITI  FIEU3  FOR  AN  ELLIPSOID 
OF  REVOLUTION  AS  A  LEVEL  SURFACE 

INTRODUCTION 

In  Volume  II  of  Stokes'  collected  papers  1/  there  are  two  mathe¬ 
matical  articles  dealing  with  attractions  and  gravity,  both  of  them 
published  in  181*9.  In  the  first  of  these,  published  in  the  Cambridge 
and  Dublin  Mathematical  Journal  and  entitled  "On  Attractions  and 
Clairaut* s  Theorem,"  Stokes  states  several  theorems*  disclaiming,  as 
he  says,  any  pretense  to  originality.  Many  of  the  theorems  are  in 
fact  due  to  Gauss. 

The  second  paper,  published  in  the  Transactions  of  the  Cambridge 
Philosophical  Society  and  entitled  "On  the  Variation  of  Gravity  at  the 
Surface  of  the  Earth,"  gives  the  formula  on  which  the  whole  gravity 
project  of  the  Laboratory  ultimately  depends.  In  the  opening  pages  of 
this  second  paper  are  two  passages  that  seem  worth  quoting  as  an  intro¬ 
duction  to  this  paper.  They  deal  with  the  general  problem  of  this 
paper  rather  than  with  the  details  of  Stokes'  Formula.  For  the  present 
purpose  the  two  concluding  sentences  of  the  first  passage  are  the  im¬ 
portant  ones. 


V  George  G.  Stokes,  Mathematical  and  Physical  Papers.  Vol.  2,  Cambridge, 

I883. 

*/  It  is  convenient  to  distinguish  between  Stokes'  Formula  given  in  this 
paper  of  181*9  and  what  is  commonly  known  as  Stokes'  Theorem,  a  very 
general  mathematical  theorem  dealing  with  the  line  integral  taken  around 
a  closed  curve  of  the  tangential  component  of  a  vector  point  function 
and  its  relation  to  a  certain  surface  integral. 
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"On  adopting  the  hypothesis  of  the  earth’s  original  fluidity, 
it  has  been  shewn  that  the  surface  ought  to  be  perpendicular  to  the 
direction  of  gravity,  that  it  ought  to  be  of  the  form  of  an  oblate 
spheroid  of  small  ellipticity,  having  its  axis  of  figure  coincident 
with  the  axis  of  rotation,  and  that  gravity  ought  to  vary  along  the 
surface  according  to  a  simple  law,  leading  to  the  numerical  relation 
between  the  ellipticity  and  the  ratio  between  polar  and  equatorial 
gravity  which  is  known  by  the  name  of  Clairaut's  Theorem.  Without 
assuming  the  earth’ s  original  fluidity,  but  merely  supposing  that  it 
consists  of  nearly  spherical  strata  of  equal  density,  and  observing 
that  its  surface  may  be  regarded  as  covered  by  a  fluid,  inasmuch  as  all 
observations  relating  to  the  earth's  figure  are  reduced  to  the  level  of 
the  sea,  Laplace  has  established  a  connexion  between  the  fona  of  the 
surface  and  the  variation  of  gravity,  which  in  the  particular  case  of 
an  oblate  spheroid  agrees  with  the  connexion  which  is  found  on  the  hy¬ 
pothesis  of  original  fluidity.  The  object  of  the  first  portion  of  this 
paper  is  to  establish  this  general  connexion  without  making  ary  hypothesis 
whatsoever  respecting  the  distribution  of  matter  in  the  interior  of  the 
earth,  but  merely  assuming  the  theory  of  universal  gravitation.  It  ap¬ 
pears  that  if  the  form  of  the  surface  be  given,  gravity  is  determined 
throughout  the  whole  surface,  except  so  far  as  regards  one  arbitrary 
constant  which  is  contained  in  its  complete  expression,  and  which  may 
be  determined  by  the  value  of  gravity  at  one  place.  Moreover  the 
attraction  of  the  earth  at  all  external  points  of  space  is  determined 
at  the  same  time;  so  that  the  earth's  attraction  on  the  moon,  including 
that  part  of  it  which  is  due  to  the  earth's  oblateness,  and  the  moments 
of  the  forces  of  the  sun  and  moon  tending  to  turn  the  earth  about  an 
equatorial  axis,  are  found  quite  independently  of  the  distribution  of 
matter  within  the  earth." 


vi 
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With  respect  to  the  determination  of  the  gravity  field  from  a 
knowledge  of  the  form  of  a  level  surface  and  the  value  of  gravity  at  a 
definite  point  on  the  surface  or  at  a  definite  external  point  Stokes 
remarks : 

"Nevertheless,  although  we  know  that  the  problem  is  always  deter¬ 
minate,  it  is  only  for  a  very  limited  number  of  forms  of  the  surface  S 
that  the  solution  has  hitherto  been  effected.  The  most  important  of 
these  forms  is  the  sphere.  When  S  has  very  nearly  one  of  these  forms 
the  problem  nay  be  solved  by  approximation." 

In  the  case  of  an  exact  ellipsoid  no  approximation  is  needed. 
Stokes  had  already  noted  the  purely  mathematical  analogy  of  the  problem 
to  a  seemingly  very  different  physical  problem,  that  of  the  steady  state 
of  temperatures  in  a  homogeneous  solid,  subject  to  certain  boundary  con¬ 
ditions.  The  solution  in  this  Part  I  is  based  on  one  of  the  commonest 
methods  of  attacking  these  problems  in  the  flow  of  heat.  The  second 
solution,  a  synthetic  one,  based  on  the  concept  of  "coating,"  is  given 
in  Part  II,  partly  as  an  illustration  of  the  convenience  of  this  concept, 
partly  because  of  certain  other  advantages.  Part  III  will  give  some 
further  developments. 
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§1  -  The  Special  Ellipsoidal  Coordinates  and  their  Properties 
The  complete  development  of  the  subject  would  occupy  an 


undue  amount  of  space.  For  the  proofs  of  many  of  the  statements  refer¬ 
ence  may  be  made  to  books  such  as  gyerly's  Fourier1  s  Series  and  Spherical 
Harmonics,  (Boston,  1895),  or  Todhunter’s  The  Functions  of  Laplace. 

Lame  and  Bessel.  (London,  1875) »  Detailed  discussions  of  the  Legendre 
functions  and  of  Laplace’s  equation  in  curvilinear  coordinate  systems 
may  also  be  found  in  Whittaker  and  Watson’s  Modern  Analysis  (Cambridge, 
1927)  or  Macfiobert’s  Spherical  Harmonics  (London,  1927). 

If  the  z-axis  be  taken  as  the  axis  of  rotat^n,  the  relation  be¬ 
tween  ordinary  rectangular  coordinates  x,  £  and  z  and  the  special  ellip¬ 
soidal  coordinates  u,  fl  and  X  is  taken  as 

x  «  c  cosec  a  Sech  0  cos  X 

y  «  c  cosec  <x  Sech  fi  sin  X  (1) 


z  «  c  cot  a  Tanh  0 . 

These  formulas  are  those  of  equation  (12)  p.  21*2  of  Byerly  with  a  re¬ 
placed  by  ^  -  a,  in  order  that  the  a  here  used  may  remain  small  in 
practice.  Byerly' s  £  and  z  are  here  interchanged. 

From  (l)  we  find 


Sech2p  ♦  Tanh*P  -  1  - 


c  cosec  a  c*cosee  a 


f 
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The  family  of  surfaces  a  -  const,  thus  represents  a  set  or  confocal 
ellipsoids  of  revolution  about  the  z-axis,  the  distance  from  center  to 
focus  being 


Vc2cosecSa  -  cScotacc  -  c. 

Equation  (2)  will  give  a  in  terms  of  x,  y  and  z.  The  equation  to  be 
solved  is  a  quadratic  in  the  square  of  a  trigonometric  function;  only 
one  root  of  the  quadratic  is  available;  the  other  root  gives  imaginary 
values  of  a.  The  result  of  solving  for  of  may  be  written 


sin' 


or 


cot  a 


i  ,  r*+c2-  y/(ra+c2)2  -  itca(xa^g) 

2(x2+y2) 

a  „  r2-c2+  \/(r^-c2)^  +  l;c*za 


(3) 


In  (3)  r2  stands  for  x2+y2+z2. 
In  a  similar  way 


cosec2u  -  cot2 a  -  1  ■ 


y* 

-It**  *  H 


c*Sech*/l  e*Sech*0  c2Tanh2p 
The  family  of  surfaces  JS  -  const,  thus  represents  a  set  of  unparted 
confocal  hyperboloids  of  revolution,  the  distance  from  center  to  focus 
being  Vc2Sech2p  +  c2Tanh3{5  -  c,  as  before.  The  ellipsoids  are  there¬ 
fore  confocal  with  the  hyperboloids  and  consequently  the  two  families 
of  surfaces  intersect  each  other  orthogonally. 

From  the  first  two  equations  of  (l) 


(W 


I  -  tan  X  (5) 

which,  for  X  -  const.,  represents  a  family  of  planes  containing  the  z-axis. 
These  planes  evidently  cut  both  the  ellipsoids  and  the  hyperboloids  ortho¬ 
gonally. 
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Any  point  in  space,  therefore,  nay  be  given  by  specifying  the 
rallies  of  a,  p,  and  X;  these  numbers  determine  the  ellipsoid,  hyper¬ 
boloid  and  plane,  the  intersection  of  which  fixes  the  point.  There  are 
apparently  four  points  of  intersection  of  the  three  surfaces,  but  the 
ambiguity  may  be  removed  by  using  the  analogy  between  the  quantities  p 
and  X  and  ordinary  latitude  and  longitude.  It  is  evident  from  (It)  that 
X  is  indeed  the  longitude  of  the  point  reckoned  from  the  xz-plane  and 
if,  instead  of  the  entire  plane  defined  by  (li),  we  use  only  the  half¬ 
plane  on  the  side  of  the  axis  specified  by  the  angle  X,  the  four  possible 
points  of  intersection  are  reduced  to  two,  one  above  and  one  below  the 
xy-plane. 

The  relation  of  the  coordinate  p  to  latitude  may  be  derived  by 
considering  the  well-known  parametric  equations  of  an  ellipse  of  semi¬ 
axes  a  and  b,  in  the  xz-plane 

x  -  a  cos  u 

(6) 

*  »  b  sin  u 

The  geometric  significance  of  the  parameter  u  and  its  relation 
to  the  vectorial  angle  Y  is  shown  in  figure  1.  The  geometric  relations 
there  indicated  suggest  the  name  "reduced  latitude"  frequently  given 
to  u,  while  by  astronomical  analogy  Y  is  the  geocentric  latitude  and  • 
is  the  geographic  latitude.  The  well-known  relations  between  these 
latitudes  are 

tan  u  ■  r  tan  Y  *  —  tan  ♦,  > 

D  a  * 

tan  Y  ■  S  tan  u  »  tan  ♦,  -  (7) 

a  aa 

tan  •  ■  r  tan  u  -  -j  tan  Y.  J 
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Figure  1  -  Meridional  Section  of  an  Ellipsoid 

a  =  semimajor  axis, 

b  =  semiminor  axis, 

PN  is  normal  to  the  ellipse,  cp  being 
the  geographic  latitude  of  P. 

The  vectorial  angle  ¥  is  the  geocentric 
latitude  of  P. 

The  angle  u  is  usually  referred  to  as 
the  reduced  latitude  of  P. 


WADC  TR  52-151 


h 


From  (2)  it  is  seen  that  the  semi-axes  of  the  ellipsoid  are 
c  cosec  a  and  c  cot  ct.  By  writing  these  for  a  and  b  in  (6)  and  intro¬ 
ducing  the  additional  coordinates  y  and  X  to  sake  the  ellipse  (6)  into 
an  ellipsoid  of  revolution  about  the  axis  we  get 

x  •  c  cosec  a  cos  u  cos  X,  > 
y  -  c  cosec  «  cos  u  sin  X,  -  (8) 

s  ■  c  cot  a  sin  e.  , 

By  comparing  (6)  and  (l)  it  is  seen  that  the  two  sets  of  equations 
would  be  identical  if  we  put 

Sech  f  *  cos  u,  'I 

[  (9) 

Tanh  fi  ■*  sin  u.  j 

The  two  equations  (9)  are  equivalent  and  mean  that  the  reduced 
latitude  u  is  the  Gudernannian  angle  of  p.  This  implies  besides  (9) 
the  additional  relations 

p  -  log*  tan  (J  ♦  |  u). 

Sink  P  »  tan  u,  *  (10) 

Cosh  p  ■  sec  m 

Equations  (9)  and  (10)  enable  us  to  use  either  u  or  jl  in  our 
equations  at  will.  It  is  found  more  convenient  to  use  p  when  obtaining 
expressions  for  the  potential  and  the  intensity  of  gravity,  but  after 
the  results  are  obtained,  it  is  sore  satisfactory  to  express  them  in 
terms  of  u.  The  fact  that  u  and  P  have  the  same  sign  and  that  u  may  be 
considered  as  the  reduced  latitude  enables  us  to  choose  between  the  two 
points  of  intersection  of  the  ellipsoid  given  by  equation  (2),  the  hyper¬ 
boloid  given  by  equation  (h)  and  the  half-plane  corresponding  to  the 
meridian  of  X  given  by  equation  (5).  The  point  on  the  positive  side  of 
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the  xy-plane  corresponds  to  a  positive  or  north  latitude  and  a  point  on 
the  negative  side  to  a  negative  or  south  latitude. 

When  x,  y  and  z  are  given,  a  is  found  from  (3);  ?  is  then  found 
by  the  following  formula  easily  deduced  from  (1) 


Sinh  ji  ■  tan  u  =  z  ■===== =L=  >  (H) 

Vx  +  y 

and  X  is  found  from  (5). 

Any  expression  for  the  potential  V  due  to  attracting  matter 
must  at  points  outside  of  such  matter  satisfy  Laplace's  equation.  For 
the  special  ellipsoidal  coordinates  a,  J3  and  X  here  adopted,  Laplace's 
equation  is 

sin2  a  +  Cosh2  fi  *  (Cosh2  {5  -  sin2  a)  ■  0.  (12) 

dor  3\ 

See  i^erly,  p.  21*2  equation  (11).  This  equation  may  also  be 
readily  obtained  from  the  equation  given  by  Whittaker  and  Watson,  p.  f>5l. 
Special  solutions  of  this  equation  are 


m-n 


n 


n 


7mn  =  (-i)  (Xjj  cos  nX  +  Bn  sin  nX)  Pm(Tanh  p)  Pm  (i  cot  a), 


n 


mn 


m+n+l  n  .  r*  d  Qm(i  cot  a) 

(i)  (An  cos  nX  +  Bn  sin  nX)  Pm(Tanh  fi)  sec  a.  ^ . —  ^ - 


In  equations  (13)  An  and  Bn  are  constants, 

Pm(Tanh  p)  represents  a  Legendre's  zonal  harmonic  of  degree  m, 

n,  .  n  d  Pm(Tanh  p) 

P„,(Tanh  pj  is  the  associated  function  Sech  3  - tt~  -  5 

d(Tanh  p) 

with  a  corresponding  meaning  for  Pm(i  cot  a). 

Qm(i  cot  a)  represents  a  Legendre's  function  of  the  second  kind. 
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The  expressions  for  Q0  (x)  and  Qx  (x),  x  being  any  variable,  are 


Q0(x)  -  \  l°ge(ll£)>  if  Ix|< 

► 

or  Q0(x)  -  \  lo2e(|§)>  if  |x|>l.  ^ 

Qx(x)  -  -1  +  |  loge^lf),  if  M<  !> 

Qx(x)  -  -i  ♦  §  i°ge(fS)>  if  lxl  >  X* 


(14) 


(15) 


From  Q0  (x)  and  (x)  the  i' s  of  higher  orders  may  be  found  by 
means  of  the  recurrence  formula 

(n+1)  Oj^x)  -  (2m+l)x  Q^x)^*  m  Q^^x)  =■  0  (16) 

The  Q® s  are  connected  with  the  P*  s  by  the  relation 


Qn(x)  -  |  Pn(x)  l°ge(j~)  "  R, 


(17) 


where 


8  ■  TT  +  WX)  *  5t£!x  pn-5<x>  + 

'  •  '  *  (&3)(n-il)  rn-2k*l(*)  +  *  *  • 


(13) 


The  series  for  R  ends  with  a  term  in  Px  (x)  when  n  is  even  and 
with  a  term  in  P0  (x)  if  n  is  odd.  If  |xj  >  1  ,  loge^il^  in  (17)  is 
to  be  replaced  by  loge 

Another  form  for  (x),  valid  when  jx|  >  1  is  the  infinite  series 


<4(*) 


ml _  1  ,  (m+1)  (m+2)  1 

( 2m+lj  (2m-l)*»*3.1  xra+1  2»(2m+3)  x®^ 

+  (m+l)(m+2)(m+3)(m+U)  _1  .  ,  # 

2.1*  (2m+3)(  2m+5)  x®^  J. 


(19) 
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From  the  definitions  we  find 


Q0(i  cot  o)  -  -  i  a, 

Qa(i  cot  a)  m  -  i  +  a  cot  a,  (20) 

Q3(i  cot  a)  -  -  ij^|  cot  a  -  (|  Cota  a  ♦  |)aj  . 

The  derivatives  of  Q0  and  Qj  will  be  needed  later;  they  are 


d  Q0(i  cot  a) 
dcf 


-  i 


(21) 


d  Qa(i  cot  «) 

-  -  i  [3  cosec  cr  (l-o  cot  ct)  -  1) 

da 


The  expressions  for  Qq  (x)  and  (i  cot  a)  satisfy  the  differential 
equations : 

(l-x2)  d3(Mx).  _  2x  +  m(m+l)  Q^(x)  -  0, 

dx2  dx 

(22) 

sina  a  - — - - — -  -  m(m+l)  Q^j(i  cot  a)  »  0* 

da3 


The  first  of  (22)  is  also  satisfied  by  the  Legendre  functions  or  zonal 
harmonics  of  the  first  kinds.  For  convenience  in  dealing  with  real 
quantities  only  let  us  define  the  quantities 

%(«)■!  QoU  cot  ®)  *  tf> 

qx  («0  ■  -  Qi(i  cot  «)  -  (l-«  cot  a),  (23) 

1 

qj  (a)  -  -  i  Q2(i  cot  a)  -  -[(3  cot2  a  ♦  1)  a  -  3  cot  a]. 


A  table  of  qjj  and  q'8  is  given  at  the  end  of  this  part,  for 
values  of  aQ  corresponding  to  values  of  the  earth's  flattening  expressed 
as  a  reciprocal.  For  the  International  Ellipsoid  of  Reference  there  is 
also  a  short  table  giving  q^cr)  and  q'-^a)  corresponding  to  certain 
elevations  above  the  ellipsoid  for  points  at  the  equator,  the  poles,  and 
at  geocentric  latitude  1*5*. 
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In  solutions  of  the  type  (13)  evidently  a  £  may  be  substituted 
for  a  £. 

§2  -  The  Expression  for  the  potential. 

The  potential  of  the  centrifugal  acceleration  U  in  rectangular  coordi¬ 
nates  is  given  by 

U  •  (x*  ♦  y3),  (25) 

ob  being  the  angular  velocity  of  rotation.  In  our  special  ellipsoidal 


coordinates  it  is 

a 

U  •  np  c3  cosec*  aSech*p,  (26) 

or  in  terms  of  Legendre's  coefficients 


2  2  2 
c  cosec  a 


[1  -  Pa(tanh  p)]. 


The  coordinate  system  is  determined  by  the  position  of  the  common 


foci.  Let  us  take  these  foci  of  the  ellipsoid  of  revolution  that  is  to 


form  an  equipotential  surface  for  the  combined  mass-attraction  and 
centrifugal  force,  and  let  Oq  be  the  value  of  a  corresponding  to  this 
particular  ellipsoid.  From  equation  (2)  it  is  seen  that  sin  aQ  equals 
the  eccentricity  of  the  meridian  ellipse,  so  that  for  the  International 
Ellipsoid  of  Reference,  for  which  the  ellipticity  is  1/297, 

Oq  -  U*  42'  IIU0I48  ■  0.0820840  radians.  Then  the  problem  consists  in 
building  up  an  expression  for  the  potential  f  »  U  ♦  ?  (V  being  the 
potential  due  to  mass-attraction),  in  which  U  and  ¥  are  expressed  in 
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terms  of  a,  fit  and  X,  which  shall  satisfy  the  following  conditions: 

(1)  V  satisfies  (12) 

(2)  V  -  o  at  an  infinite  distance  from  the  origin, 

(3)  W  ■  (U  ♦  V)  is  independent  of  and  X,  when  a  is  put  equal 
to  Oq}  that  is,  since  V  is  a  function  of  a0  only,  W  is  constant  on  the 
ellipsoid  a  *  o0t  which  is,  therefore*  an  equipotential  surface. 

(ii)  V  must  approximate  to  ~  sis  r  increases,  k  being  the  gravi¬ 
tation  constant,  M  the  mass  and  r  the  radius  vector. 

Condition  (l)  is  satisfied  by  any  solution  formed  of  special 

solutions  of  either  type  shown  by  equation  (11).  The  semi-axes  of  the 

coordinate  ellipsoid  being  c  cosec  a  and  c  cot  a,  it  is  seen  that  a  •  o 

corresponds  to  a  point  at  infinity,  and  that  only  solutions  of  the 

second  type  under  (13)  may  be  used  if  condition  (2)  is  to  be  satisfied, 

for  the  Pm(i  cot  a)  of  the  first  type  increases  without  limit  when  a 

approaches  zero,  whereas  as  may  be  seen  from  the  series  (19)  by  putting 

i  cot  a  for  x,- — c approaches  zero  with  a.  It  is  clear  that 
d  (i  cot  <0 

since  the  equipotential  surface  is  to  be  symmetrical  about  the  axis  of 
rotation,  there  can  be  no  terms  involving  the  longitude,  that  is, 
n  »  o  in  all  cases,  and  since  the  surface  is  symmetrical  about  the 
equatorial  plane,  there  must  be  no  terms  that  change  sign  with  a* 
this  condition  excludes  from  V  terms  in  which  m  is  odd.  Sty  using  (27) 
and  noting  that  PQ  (Tanh  fi)  ■  1  we  therefore  can  write 
w  -  V  *  V  -  <-c-23c-03eC-  a  ♦  Ao0o(i  cot  a) 

♦  Pa  (Tanh  p)  [A,  0,(1  cot  a)  -  “f-Ji^cosec*  a]+A4P4(Tanh  j*)Q4(i  cot  a)  (28) 
+  A£P6(Tanh  0)  Q^(i  cot  «)♦•••. 
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The  powers  of  i  necessary  to  make  the  terms  real  are  supposed 
to  be  included  in  the  A' s.  If,  according  to  condition  (3),  W  is  to 
be  independent  of  0  when  a  -  Oq  then  o  •  A4  ■  A6  -  As  -  •  •  •  . 


A^U  cot  Oq)  - 


2  2  2 
or  c  cosec  a. 


2-o 


which  determines  Aa.  Accordingly  (28)  becomes,  after  recombining  the 
separate  parts  of  the  expression  for  the  centrifugal  force  into  their 
original  form  (26),  and  substituting  for  Qa  in  terms  of  qa 


W  -  cot  o)  ♦ 


2  2  2 
aa  c  cosec  a 


Sech3  0 


2  2  2 
or  c  cosec  a„ 


Pa  (Tanh  fl)  qa(a). 


We  can  satisfy  condition  (4)  by  choosing  A0  suitably;  to  do  this 
we  must  ascertain  the  relation  between  a  and  r,  when  r  is  large  or  a 
small.  By  squaring  and  adding  equations  (1)  we  get 

r  -  cVcosec2  o  Sech3  fl  +  cot2  a  Tanh3  .  (3] 

But  the  expansions  for  cosec  a  and  cot  *  in  power  series  are 

la  la 

cot  a  -  —  —  rr  •  *  •  and  cosec  a  -  —  *  t  '  *  •• 
a  j  a  o 

These  substituted  in  (31)  give 


r  *  °V  j  Tanh3  *  j  Sech2  j!  ♦  terms  in  a3,  etc.  , 

from  which  it  appears  that  r  -  —  is  a  good  approximation  when  a  is  small 
or  r  large. 

By  equation  (1?)  Qa(i  cot  a)  or  q^a)  is  of  the  third  order  when  a  is 

small,  a  being  of  the  first,  so  the  term  in  (30)  which  involves  qa(a) 

need  not  be  considered  in  satisfying  condition  (4).  It  is  thus  seen  that 

R  M  i 

condition  (k)  is  satisfied  by  taking  A0  - - —  and  the  expression  (28) 
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for  W  becomes 


k  M  a  u?  c2  cosec2  a 

W  -  - +  - 


Sech2  fi  * 


2  2  2 
a>  c  cosec  «o 


3  q,(«o) 


Pa(Tanh  p)  q2(ct).  (32) 


This  equation  satisfies  all  conditions  and  is  therefore  the 
potential  function  that  has  the  ellipsoid  of  revolution  ec  •  aQ  as  an 
equlpotential  surface  for  the  combined  effects  of  mass-attraction  and 
centrifugal  acceleration. 

§3  -  The  formula  for  the  Intensity  of  Gravity. 

As  a  preparation  for  obtaining  the  gravity  formula  let  u3  derive 
the  distance  between  the  two  ellipsoids  characterized  by  the  coordinates 
a  and  a  ♦  Aa,  A«  being  infinitesimal.  Take  a  point  P  on  the  first  ellip¬ 
soid  whose  coordinates  are  ( a,  p,  X),  or  in  rectangular  coordinates 
(x,  y,  z),  and  a  point  Q  on  the  second  whose  coordinates  are  (a+  Act,  p,  X), 
or  in  rectangular  coordinates  (x  ♦  Ax,  y  +  Ay,  z  +  Az).  Since  the  system 
of  coordinates  is  orthogonal,  the  infinitesimal  line  of  length 

P Q  *  VAx2  +  Ay2  +  Az2  will  be  perpendicular  to  both  ellipsoids  and  will 
be  the  distance  required.  Since  a  alone  varies 

Ax=*  A«  +  terms  of  higher  order. 


Ay-  |i«  ♦  terms  of  higher  order, 

Az  -  ~  Act  +  terms  of  higher  order. 

Or  by  differentiating  the  equations  (1),  squaring,  adding  and  simplifying 

PQ  *  \/Ax*  +  Ay5*  +  Az2  -  c  cosec  a  Vcot2  a  +  Tanh2  p  A«  ♦  terms  of 
higher  order, or  since  PQ  is  an  element  of  the  normal  n, 

An  *  -|£  Aa  »  c  cosec  a  Vcot2  a  +  Tanh2  p  Aa  +  terms  of  higher  (33) 

<7\X  * 

order.  Since  a  decreases  outward  from  the  origin,  the  above  expression 
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implies  that  the  normal  n  is  directed  inward. 
For  brevily  let  us  put 


h  ■  c  cosec  a  vcot2  a  ♦  Tanh2  £  ,  (3l 

so  that  equation  (33)  is  equivalent  to 

3n  , 

3^  "  h 

It  will  be  convenient  to  have  the  special  value  at  the  ellipsoid 
a  -  aQ.  If  we  call  a  the  semi-axis  of  the  ellipsoid, 
a  -  c  cosec  Oq  and  e  -  sin  a0.  Therefore  by  (8)  and  (9),  equation  (32) 


gives,  after  a  slight  reduction 


tan  aQ 


ho  a  VI  ♦  ban*  aQ  sin2  u 

The  subscript  o  attached  to  h,  indicates  the  value  for  the 
ellipsoid  a  •  o0. 

If  we  introduce  the  geographic  latitude  #  instead  of  the  reduced 
latitude  u,  we  find  for  ho  after  some  transformation  by  the  help  of 
equations  (7) 

h„  -  .Sa...—  .  (36! 

VI  -  sin  a0  sin  ♦ 

The  value  of  the  intensity  of  gravity  (g)  at  any  point  may  be 
found  by  dividing  the  difference  in  potential  between  the  two  ellipsoids 
at  an  infinitesimal  distance  from  one  another  by  that  distance*  and 
dropping  the  terms  of  higher  order  which  vanish  when  the  limit  is  taken. 

1.3W 
6  *  I  jff 

(37) 

m  H  4 
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By  differentiating  (32)  to  get  3—,  and  using  (2it)  we  find 


kM  s 
"[c"  “  0 


c3  cosecaa  cot  a  Sech3p  +  —  *a(Tanh  P) 


_c\/cot3a  + 


sin  a 


TanhzjJ 


This  is  the  expression  for  the  intensity  of  gravity  at  any  point 
(a,  p,  X)  outside  the  ellipsoid  a  =  a0.  The  case  when  the  point  is  on 
this  ellipsoid  is  included  as  a  limiting  case,  and  for  this  case  (38) 
can  be  simplified. 

Let  us  put  o  *  «o  and  introduce  the  abbreviation  F  defined  by 

„  <&(ao)  sin  ao 

F '  3  (3S 

In  the  first  expression  F  is  indeterminate  in  form  when  a*0,  and  this 
expression  is  therefore  inconvenient  for  purposes  of  computation  when 
aQ  is  small,  as  in  the  case  of  an  ellipsoid  approximating  the  figure  of 
the  earth.  In  this  case  it  is  convenient  to  use  the  series  (19)  for 


Q2(i  cot  a),  which  gives 


q2(a)  *  -i<32(i  cot  a)  mjr  [ tan3a-  y  tan5a+^  tan7o-  ~  tan9a  •  •  •] 
*  — ^  tan3a-  ^^tansa  +  tan7 a  -  tan9 a  •  •  • 


qj^(a)  tan2 a  -  29  tan4a  +  £  tan6<t  -  ^  tan6a  *  •  •], 

tan2 a  I4.  tan4a  ,  6  tan6a  8  tan8a 

“  sec  ccl - 5= - j - +  — - -£j—  .  .  .] 

so  that  10  .  5  .  20  , 

sec  a0[l-  7  tan  cCq+3  tan  oQ-  tan°a0  +  •  •  •] 

F  » ■■■■ - 1  r - “S - - 55 - 7 -  •> 

1-  y  tan2a0  +  2.  tan4a0  -  yy  tan6o0  +  •  •  * 

or  approximately  by  division 

F  -  sec  aQ[l-  2  tan2a0+  Jji  tan4a0-  J±|i|  tan^ac+  •  •  •]• 
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Putting  a  »  aQ  in  (38),  introducing  F,  substituting  for  p  in  terms  of 
u  and  simplifying  gives  for  g0,  the  value  of  gravity  on  the  ellipsoid, 

g0  ■•F-y  -  6>aa  cos  o0  cos2  u  ♦  &>2a  F  Pa(sin  u)  t  4  n.o\ 

'  a  *  Vl  -sin2  «*0  cos2  n 

Instead  of  the  constants  M  and  u  let  us  express  the  result  in 

terms  of  ge,  the  value  of  gravity  at  the  equator,  and  m  the  ratio  of 

the  centrifugal  acceleration  at  the  equator  to  gravity  there,  or 

m  -  w2  a/ge. 

By  putting  u  -  o,  and  w2  a  «  i  g0  in  (1*2)  and  solving  for 
k  M/a2  we  find 

-  ge  [cos  aQ  +  in  (i  F  cos  oQ)].  (1*3) 

d>  ^ 

By  using  this  (1*2)  reduces  to 


g0  *  [1  *  m  (2  F  sec  <xQ  +  1)  sin2  u]«  (1*1*) 

Vl  +  tan  a0  sin2  u  £ 

This  is  an  exact  expression  for  gravity  on  the  surface  of  the 
ellipsoid  at  a  point,  the  reduced  latitude  of  which  is  u,  sin  ccQ  being 
equal  to  e,  the  eccentricity  of  the  ellipsoid  and  F  being  given  by 
(39),  (1*0)  or  (1*1).  This  can  be  expressed  directly  in  terms  of  the 
geographic  latitude  $.  From  (7) 


sin  u 


sin  $  cos  <r„ 


Vl  -  sin2  a,-,  sin2  * 
and  after  substituting  this  in  (1*1*)  and  reducing  we  get 

(^  F  +  cos  %)  sin2  *  ] 


So  m  Se 


Vl-  sin2  an  sin2  # 


1  +  m  cos  a„ 


1  -  sin2  Oq  sin2  # 
This  expression  is  exact  for  an  ellipsoid  of  any  eccentricity. 


XUS) 
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On  the  assumptions  stated,  the  expression  for  g0  is  rigorous. 
Moreover  it  is  in  finite  terms,  for  the  quantity  cos  a0(^  F  +  cos  Oq) 
is  in  any  given  case  merely  a  numerical  constant,  although  the  way  in 
which  this  quantity  involves  o0  is  somewhat  complicated. 


§U  -  Various  Deductions  from  the  Formula  for  Gravity. 
Clairaut1  s  Equation. 

By  putting  $  *  90*  in  (1*5)  or  u  -  90*  in  {kh)  we  get  for  gp, 
the  value  of  gravity  at  the  pole 


gp  ■  ge  cos  «o  j 


1  +  m(|  F  sec  a0  +  1)1  , 


and  by  forming  the  expression  for  the  quantity 


m 


we  find 


a  -  hCSi 

P  8e 


P 


-  {1  f  +  cos  a0)  m  -  (1  -  cos 
ge 


ao) 


(1*7) 


Let  us  introduce  the  flattening  f,  which  is  connected  with  Oq  by  the 
relation 

e2  -  sin2  aQ  *  2f  -  f2.  (1*8) 

Let  us  consider  m,  e2  and  f  as  small  quantities  of  the  same  order 
of  magnitude;  as  a  matter  of  fact  m,  being  equal  to  about  1/288,  is  very 
nearly  equal  to  f.  If  in  ( i*7)  we  consider  only  quantities  of  the  first 
order,  we  get  by  using  (1*1) 

P  -  -  |  m  -  f .  (]*9) 

This  is  the  familiar  form  of  Clairaut* s  equation  and,  though 
accurate  to  the  terms  of  the  first  order  only,  is  sufficient  for  many 
purposes.  It  is  desirable,  however,  to  include  terms  of  higher  order. 


•151 
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By  putting  in  (1*1)  sec  Oq  »  y/l  *  tan2  cc  »  1  ♦  ~  tan*  «o  -  ^  tan4  •  •  • 
multiplying  out  this  series  and  the  other  series  in  (1*1)  and  substituting 

*“*  ■  fl§F  -  *  f  ♦  3  f  ♦  k  f  ■  •  •  , 

also  noting  that  1  -  cos  Oq  *  f,  we  find  after  reduction 

f»  *  «§m(l-^f-^  f*— )-f-  §a*(f)-f  (SO) 

the  function  x(f)  being  defined  by  the  series  in  parentheses. 

This  is  Clairaut's  equation  extended  to  terms  of  even  higher 
order  than  are  needed  in  practice. 

Equation  (f?0)  applies,  of  course*  only  when  the  level  surface  is 
an  exact  ellipsoid  of  revolution. 

Relation  Between  the  Newtonian  Constant  k  and  the  Mean  Density 
of  the  K&rth. 

The  mean  density  Pm  is  defined  by  the  condition  that  the  mass  of 
a  homogeneous  ellipsoid  of  the  adopted  dimensions  shall  have  the  same 
mass,  that  is,  that 

k  -rra2  b  Pm  -  ^  -n-a3  cos  aQ  Pm  «  M  (f?l) 

Substituting  this  value  of  M  in  (1*3)  and  expressing  F  and  aQ  in  terms 
of  the  flattening  f  we  get 

k/Ja  «  [1  ♦  \  ■  (1  ♦  I  f  ♦  rr|  f2  •  *  •)).  (S2) 

*t  wa  c  I  hai 

An  equivalent  form  is  given  without  proof  by  Cassinis  in  the 
Bulletin  G^odesique,  No.  26,  April-May-June,  1930,  p.  1*0.  Cassinis 
gives  also  an  additional  term  for  *(f)  following 

-  2^  f2,  namely  -  f3.  The  correctness  of  this  additional 

term  has  been  verified  by  the  writer. 
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The  Mean  Value  of  Gravity 

It  is  understood  that  the  mean  is  to  be  taken  with  regard  to  area. 
The  value  of  mean  gravity,  is  therefore  given  by 

kas 

g"‘  5^ 

Where  dS  is  an  element  of  area  of  the  ellipsoid  and  the  integration 
covers  the  entire  area. 

The  reduced  latitude  u  is  convenient  in  evaluating  these  integrals. 
In  this  case  we  have  dS  *  2  -Trxds  where  ds  is  an  element  of  meridional 
arc.  From  x  ■  a  cos  u,  y  »  b  sin  u 


we  find 


dS  «  2  'Jr ab  cos  u  vl  +  tan2  aQ  sin2  u  du. 


As  is  well  known,  this  may  be  evaluated  in  closed  form  in  terms  of  the 
elementary  functions  but  the  expression  obtained  is  inconvenient  to 
calculate  when  tan  a0  is  small.  It  is,  therefore,  simplest  for  the 
purposes  of  practical  computation  to  expand  v4  +  tan*  «cQ  sin2  u  in 
powers  of  sin  u  and  integrate  term  by  term. 

In  this  way  we  find 

S  -  T  dS  «  li  nr  ab  [1  +  ^  tan2  a0  -  tan4  aQ  +  tan6  (53) 

J  -tt 
u«__ 

In  calculating  ^  g0  dS  it  is  convenient  again  to  expand 

v4  +  tan2  aQ  sin2  u  in  powers  of  sin  u.  In  this  way  by  using  equation 

(liii)  for  gQ  we  find 
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C  \0dS  -  2-rrabge  f^[l+m(|  F  sec  Oq  ♦  1)  sina  u]  cos  u  du 
J-rr  ^  Tf 


■  li«abge  [1  +  j  F  sec  Oq  +  1)] 

-  Wbge  [1  ♦  |  m(l+  ^  fcanaaQ  -  ^  tan4a0  •  .  •)]. 


Therefore 


-  jgo*33  ,  ^4  mU  ♦  ^  tanaa0  -  tan4^  *  ♦  y 

S  1  +  £  tan2Oo-  jj^tan4a0  ♦  tan6ot0  •  •  • 

-  g0  C1  “  Z  ^““o  4  ^7?  tan4<Io  -  tSto  tan6a0  •  •  . 


^  tan  °o  -  jgZo  tan*aQ 


|“(1+  ^>tan2co  -  ifsfe  t*n4ac 


ll-if- 


The  value  gm  of  mean  gravity  is  of  interest  because  in  practice 
it  is  usually  better  determined  than  gravity  at  the  equator.  Formulas 
based  on  rather  different  values  of  gravity  at  the  equator  and  on 
rather  different  flattenings  often  give  nearly  the  same  value  for 
mean  gravity. 
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§  5-  Other  Forms  of  the  Expression  for  Theoretical  Gravit 


The  quantity  ~  F  cos  aQ  +  cos*  «o  depends  merely  on  the  flatten- 
ing  of  the  ellipsoid.  For  small  flattenings  its  value  is  approximately 
For  brevity  let  us  denote  this  quantity  ^  F  cos  Oq  ♦  cos*  a  by  C. 

The  value  of  C  is  obtained  from  (lj.0)  by  multiplying  by  cos  a  and 
dividing  the  denominator  into  the  numerator.  We  find  in  terms  of  tan  a. 

F  cos  oQ  «  1  -  |  tan*  Oq  +  tan4  a 0  -  tan6  <*©•••  (56) 


Since 


_ _ 2  _  sin*  ao 

tan  an  *  1  ■■  ■  ■■ 

°  1-e  in2  ar 


«  e*  +  e4+  e*  •  •  • 


We  find  by  substitution 


F  cos  a  *  l  -  k  e*  -  e4 - 9^  •  •  • 

0  7  Wf  11319 

We  have  in  terms  of  the  flattening  f  by  means  of  the  relation 


e2  -  2  f  -  f2. 


F  cos  aQ  -  1  -  7  f  +  ^  f2  +  u.319  f3  *  *  • 


The  forms  ($8)  and  (60)  are  more  convergent  than  (56). 
In  terms  of  e2  and  f  we  have 

C«!-i2e2-4"e4-  e6  .  .  . 

2  7  3773 

5  26  .  ^  59  „2  ^  _32_  - 

2  "  7  f  C9f  3773  ^  *  * 

Returning  to  (U5)  we  write  it 


g0  -  ge  -1  Vl-e2  sin5*  #  + 


m  C  sin  $ 
V'l-e^  sin2  < 
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This  form  suggests  an  expansion  in  powers  of  sin2#;  we,  therefore, 

assume 

So  “  Set1  +  C2  sin2*  +  c4  sin4#  +  C6  sin6*  *  •  •].  (64) 

By  expanding  the  radicals  in  (63)  and  noting  that  m  is  a  quan¬ 
tity  of  the  same  order  as  e2  we  find 


The  numerical  coefficient  of  C2k  is  that  of  jr  in  the  expansion 
of  vl  -  y.  The  numerical  values  of  the  C's  will  be  given  in  §6. 

For  actual  computation  of  tables  the  expansion  in  multiple  angles 
is  probably  as  convenient  as  any.  If  we  put 

go  *  !  B0  ♦  B2  cos  2#  +  B4  cos  4#  +  B6  cos  6#  *  •  •  (67) 

and  use  the  expressions  for  even  powers  of  the  sine  in  terms  of  the 
cosines  of  multiple  angles,  namely, 

sin2#  ■  cos  2# 

sin4#  »  ^  cos  2#  +  ^  cos  4# 

sin6#  -  ^  ^  cos  2#  ♦  cos  4#  “  -yg~  cos  6# 

sin8#  -  ^  ^  cos  2#  ♦  -^2  cos  cos  *  +  125  cos  8# 
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we  get  by  collecting  coefficients 


\  Bo  "  ge^1+  \  cz  +  |  c4  +  iEC6  + 

B2  *  "Se(|c2  +  |C4+  llC6+  ftc6  *  *  •) 

(68) 

B4  -  ge(|c4  +  ^C6  +  ^C8*  •  •) 

B6  *  “Se(j2  c6 +  5^  c8  *  *  * ) 

B8  3  ««<I2&V  *  •) 

The  coefficients  of  the  expansions  of  (l-e2sin2*)  c  may  be 
expressed  as  hypergeometric  series  in  terms  of  e2  or 

1  -  vl  -  e2  2  ao 

n *  i+  =  tan  *2“ ,  which  latter  series  are  even  more  rapidly  con¬ 

vergent  than  those  in  e2.  The  B«s  of  equations  (68)  may  also  be  ex¬ 
pressed  as  elliptic  integrals.  The  functions  qx  and  q2  are  also 
expressible  as  hypergeometric  series,  but  for  the  purpose  in  hand  the 
elementary  treatment  here  given  seems  adequate. 

Cassinis  ^/puts  the  formula  for  gravity  in  a  form  resembling 
that  previously  in  use,  namely 

gQ  *  ge(l+P  sin2*-px  sin22*-j32  sin2*  sin2 2* 

-jJ3  sin4*  sin22«  •  •  •), 

a  form  in  which  all  terms  after  the  term  in  0  contain  the  factor 
sin2 2*  and  thus  vanish  both  at  equator  and  pole.  It  will  be  seen  that 
this  3  is  the  same  as  the  £  of  equations  (hi)  and  (5>0) .  Helmert' s 
formula  of  1901  stops  with  the  term  in  (3*  5 /.  As  a  matter  of  fact  in 

G."  Cassinis'.  Bulletin  G6od6sique,  No.  26,  April-May-June,  1930,  p.  hO 
5 /  Ber.  d.  Kbn.  Pr.  Akad.  d.  tfiss.,  XIV  1901,  S.328, 
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the  present  instance  the  succeeding  terms  are  very  small.  The  values 
of  the  {Ms  can  be  determined  in  terms  of  the  C's  previously  used  by 
means  of  the  relation  (68). 

By  substituting  this  in  equation  (69),  collecting  terms  and 
comparing  coefficients  we  find 


p  «  ca  +  c4  +  c6  +  ca  •  *  *  . 
J-fcPx-Ca. 


(70) 


h  -  J(c4  +c6  +  ca  *  *  *)• 

Pi  “  p2 )  ■  , 


(71) 


h  =  J(C6  +c8  *  *  *)• 

WPa-P3)  “  C6» 


(72) 


or 

P3- j(C8  •  •  •)  (73) 

Another  form  given  by  Cassinis  (l.a.c) 

By  reducing  to  a  common  denominator  within  the  braces  and  using, 
as  before,  the  abbreviation  C  for  |  F  cos  aQ  +  cos2aQ)  (1±5)  may  be 
written 


So 


ge[l.t  (mC  ~  e2)sin2S] 
Vl-e2  sin2* 


(7k) 


It  can  be  shown  that 


mC  -  e2  =  p-f-fp, 

so  that  (7k)  may  be  written  in  Cassinis* s  form 

get1 + (P-f-fp)sin2*]  ge[l+ (P-f-fp)sin2*] 

m  -  -  _  -  - - - - - - -  m  _  -  -  -  . 

\&-Vsin-*  Vl  -  f  ( 2  -  f)sin-* 


(75) 

(76) 
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The  quantity  W  =  1  -  e2sin2*  is  given  in  the  tables  for  the  Inter¬ 
national  Ellipsoid  of  .Reference  &]  from  **0*  to  **U5°  and  of  the 

quantity  V  =  ^  '$  *  Vl  -  e^  W  from  *  *  h5“  to  =  90* .  But  if  these  tables 
<* 

are  to  be  used  for  the  systematic  calculation  of  g0,  it  is  simpler  to 
use  a  transformation  involving  the  "great  normal"  1/  N,  a  transformation 
that  avoids  the  discontinuity  in  tabulation  at  4*  and  also  the 
necessity  of  passing  from  W2  to  •».  By  using  the  relation  between  N  and 
V  or  W  we  may  write 


go-ge|[l+  (p-f-fP)sin2$] 
-gef  [1+  (nC  -  e2)sin2*]. 


(77) 


The  quantity  gg/a  is  constant  and  the  logarithm  of  N  may  be  taken 
directly  from  the  tables. 

A  symmetrical  formula  involving  the  shape  of  the  ellipsoid  and 
the  values  of  gravity  at  Equator  and  pole,  or  ge  and  gp,  respectively, 
has  been  given  by  Somigliana  &/ .  Let  denote  the  ratio  of  the  semi- 
major  to  the  semi-minor  axis  or  a/b  *  l/\A-  e2  then  Somigliana  finds 

g2( jZc°s2* +  sin2*)  -  (j  ge  cos2*  +gp  sin2*)2.  (73) 

A  similar  formula  may  be  written  in  terms  of  the  reduced  latitude  u. 

It  reads 

g2(  j2sin2u+ cos2u)  -  (gecos2u  + j  gp  sin2u)2.  (79) 


bj  Tables  de  1*  Ellipsoide  de  Reference  international - calculees  sous 

la  direction  du  &<£ndral  G.  Perrier - par  E.  Hasse,  Paris,  1928. 

7/  This  quantity  is  also  the  radius  of  curvature  in  the  prime  vertical. 

&/  C.  Somigliana.  Atti  della  fteale  Accademia  Nazionale  dei  Lincei, 

Anno  CCCXXIV  (1927),  Series  6,  Vol.  £,  p.  319. 


2k 
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The  proof  of  the  first  may  be  taken  from  Somigliana' s  paper*  or  the 
proofs  of  both  may  readily  be  derived  by  the  reader. 

§6  -  Numerical  Values 

For  the  International  Ellipsoid  the  value  of  the  semi-major  axis 

a  is  6  378  388  meters  and  the  value  of  the  flattening  f  is  1/297.  From 

these  we  find 

log  a  -  6. 80U  71093 
log  f  -  7.527  21^355  -  10 

The  period  T  of  rotation  of  the  earth  from  star  to  star  is  86164.09890 

mean  solar  seconds.  This  period  is  nearly  the  same  as  the  "sidereal*' 

day  so-called  expressed  in  mean  solar  seconds  but  should  be  used  rather 
than  the  so-called  sidereal  day  in  computing  the  centrifugal  acceleration. 
The  so-called  sidereal  day  is  really  the  equinoctial  day  or  interval 
between  successive  transits  of  the  vernal  equinox.  A  more  appropriate 
name  for  the  so-called  sidereal  day  would  be  the  equinoctial  day.  From 
this  value  we  find 

u  -  2t/T  -  0.000  07292  11515 

log  w  -  9.8628  53518-  10 

The  Association  of  Geodesy  on  the  recommendation  of  its  Gravity 
Committee  adopted  as  the  value  of  ge,  or  gravity  at  the  fcquator  at 
sea  level,  978.049  cm/sec8.  This  value  is  based  on  the  work  of  Heiskanen  sj. 
With  these  values  of  a,  «,  and  ge  we  find  for  m  -w3  »/ge 

log  m  -  7.5400  57356  -  10 

9/  W.  Heiskanen.  1st  die  Erde  ein  dreiachsiges  Ellipsoid?  Ger lands 
Beitrage  zur  Geophysik,  Vol.  X1a.  1928,  p.  jpb. 


>7 ADC  TR  52-151 


25 


With  these  values  as  a  basis  we  find  for  the  various  quantities  for 
which  formulas  have  been  given 


C  -  2.1*87  50763  8 
Ca  -  0.005  261*90  98 
c4  -  0.000  02332*  61* 

C6  -  0.000  00012  72 

Cg  -  0.000  00000  07 

P  •  0.00 5  28838  1*1 
px  *  0.000  00586  86 
p2  -  0.000  00003  20 

x(f)  -  0.998  361*55  21 
1 

2  B0  -  980.632  272  cm/sec2 
B2  *  -  2.586  11*5 

B4  -  +  0.002  878 

Bfc  «  -  0.000  001* 

log(p-f-fp)  -  7.279  5699  -  10 
log  ge/a  "  6.185  61*96  8  -  10 

We  have,  therefore,  the  following  working  formulas  for  computing  the 
acceleration  of  gravity  at  sea  level  in  latitude  *,  the  unit  being  the 
gal  (cn/see2) 

g0  -  978. 0l*9[  1+0.005  261*91  sin2* +  0.000  02335  sin4* 

♦  0.000  00013  sin6*]. 

gQ  -  978. 0l*9[  1  +  0.005  28838  sin2* -0.000  00587  sin2 2* 

+  0.000  00003  sin2*  sin2 2*]. 

gQ  -  980.632  272  -  2.98611*5  cos  2* 

+  0.002  878  cos4*  -  0.000  001*  cos6*. 

g0  -  [6.185  61*968  -  10]  N  -fl  +  [7*  2795699  -  10]  sin2*}  . 


(A) 

(B) 

(C) 

(D) 
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In  conformity  with  a  usual  convention  among  astronomers,  the 
square  brackets  in  the  formula  (D)  indicate  that  the  number  within  is 
the  common  logarithm  of  the  coefficient  concerned. 

§7  -  The  q-functions.  Tables  and  Computation 

Expression  (32)  for  the  potential  and  expression  (38)  for  gravity 
at  any  point  in  space  external  to  the  ellipsoid — and  including  the 
surface  of  the  ellipsoid  itself — are  exact.  At  points  sufficiently  far 
from  the  surface  of  the  ellipsoid,  the  concept  of  a  field  of  force 
representing  the  combined  effects  of  Newtonian  attraction  and  centrifugal 
force  becomes  meaningless,  but  in  practice  the  points  likely  to  be 
considered  are  far  within  the  limit  of  applicability. 

These  equations,  (32)  and  (38),  were  derived  with  the  aid  of 
curvilinear  coordinates  based  on  three  families  of  mutually  orthogonal 
surfaces,  one  family  being  a  set  of  ellipsoids  of  rotation.  But  of  all 
these  ellipsoids  the  only  one  that  is  of  special  interest  to  us  is  the 
ellipsoid  a  ■  aQ,  where  sin2a0  -  ea  -  square  of  eccentricity  of  meridian 
ellipse.  Tne  condition  a  <  corresponds  to  ellipsoids  of  the  same 
family  outside  the  ellipsoid  of  reference  for  which  a  •  aQ.  Unfortunately 
these  outer  ellipsoids  do  not  give  us  much  help  in  formulating  our  ideas 
about  the  level  surfaces  outside  of  a  -  aQ.  The  ellipsoids  of  increasing 
axes,  corresponding  to  decreasing  values  of  oc,  become  less  and  less 
flattened  as  we  recede  from  the  center.  On  the  contrary,  the  actual 
level  surfaces  become  more  and  more  flattened  under  the  same  conditions.10/ 


l0/  This  is  due  to  the  inclusion  of  the  centrifugal  force.  The  level 
surfaces  of  the  Newtonian  attraction  alone  would  become  more  and  more 
nearly  spherical. 
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Moreover,  the  level  surfaces  are  not  exact  ellipsoids  but  spheroids 
depressed  in  Diddle  latitudes  below  the  ellipsoids  having  the  same 
equatorial  and  polar  axes. 

It  would  seem,  therefore,  that  expressions  (32)  and  (38)  would 
seldom  be  used,  so  that  the  values  of  q2(a)  and  q’(a),  where  a  *=  aQ  t 
could  be  computed  as  occasion  might  require.  Thus,  it  would  also  seem 
that  elaborate  tables  of  these  functions  would  be  of  little  use. 

But  in  order  to  give  an  idea  of  the  general  trend  of  these 
functions,  two  tables  have  been  prepared.  Table  I  gives  the  values  of 
q2(a0)and  ^(a^for  values  of  aQ  corresponding  to  ellipticities 
(flattenings)  denoted  by 

f  =,  2  sin2  -  — - —  — - — _ i 

2  296.0  »  296.1  *  293.5  * 

This  range  covers  all  modern  ellipsoids  of  reference  likely  to  be  of 
interest.  For  obvious  reasons,  figures  for  three  special  values  outside 
this  range  are  given,  namely  for  the  ellipsoids  of  Clarke  (1880), 

Clarke  (1866),  and  Bessel. 

Table  I  applies  to  the  ellipsoid  of  reference  only,  as  is 
indicated  by  the  subscript  zero. 

Table  II  is  for  the  International  Ellipsoid  only  (flattening 
equals  I/297  exactly)  but  applies  to  a  region  above  (outside)  this 
ellipsoid.  The  argument  is  the  altitude  of  a  point  above  the  equator. 
The  table  gives  the  values  of  a,  q*(a)  ,  and  q2(a)  for  radii  exceeding 
the  radii  of  the  earth  by  10,  20,  30,  •  •  •  100  kilometers  and  at  three 
points,  the  equator,  the  pole,  and  geocentric  latitude  1*5“ .  At  the 
equator  and  pole  the  excess  radii,  10,  20,  30  •  •  •  kilometers,  are 
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altitudes,  but  clearly  at  latitude  U5*  the  excess  radius  is  not  exactly 

an  altitude,  whatever  definition  we  may  give  to  the  concept  altitude. 

If  Ar  is  the  excess  radius,  the  difference  between  Ar  and 

f* 

altitude  defined  in  ary  reasonable  way  is  of  the  order  of  —  Ar  at 
latitude  45*,  the  latitude  at  which  ambiguities  in  the  definition  of 
altitude  produce  the  largest  discrepancies.  Possible  definitions  of 
altitude  will  be  considered  in  Part  III. 

The  second  of  equations  (3)  is  readily  transformed  into 


tan  a  « 


[r  -  cJ)*  ♦  ire4**  ♦  ra  -  c* 


r  -  c'r  ♦  kc  z  - 


(r2-  ca) 


For  the  equator  we  may  put 


and  for  the  pole 


whence  we  find: 
at  the  equator 


x  «  r,  y  «  z  «  0, 


z  -  r,  x  -  y  -  0, 


sin2  a  m  — 


tan2  a 


2  2  * 

r  —  c 


at  the  pole 


sin2  , 
r  +  c 


tan2  a  -  . 

r2 
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With  the  usual  notation  for  the  semi-axes  of  the  ellipsoid  of  reference, 
c2  -  a2e2.  Equations  (80)  and  (8l)  give  for  the  surface  of  the  ellipsoid 
at  ooth  equator  and  pole 


sin2ar 


» 


tan2a0  *  c2 


y 


wnich  are  merely  our  adopted  definitions,  as  they  should  be. 

To  determine  a  for  a  point  outside  the  ellipsoid  of  reference, 
it  may  often  be  convenient  to  expand  the  radicals  in  (3)  or  (3a)  in  series 
after  taking  r2  ♦  c2  or  r2  -  c2  outside  the  radical.  Since  all  meridian 
ellipses  are  alike,  we  put 


7-0 

x  *  r  cos  Y  , 
z  =>  r  sin  f  , 

where  Y  is  the  geocentric  latitude  in  the  expression  for  sinaa. 
For  brevity  call 


sin2 


"  r2  ♦  c2 


(82) 


Then  from  (3)  we  find 

sin2«  »  sin2 Op  +  sin4ap  cos2Op  cos2f  +  2  sin^Op  cos4Op  cos4Y 

(83) 

♦  5  sin8Op  cos6<ip  cos6!  +  li;  sinl0Op  cos8Op  cos*Y  +  •  •  •. 


From  (3a)  we  find,  putting 


tanaa. 


e  a  2  * 
r  —  c 


t*n*a  •  tan*ae  -  tan4ag  secaa0  sin2Y  ♦  2  tan6^  sec4ag  sin4Y 
-  5  tan8ae  sec^sin6!  +  lij  tanl0ae  sec8ae  sin8Y  +  •  • 


(810 

(85) 
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The  coefficients  2,  ±  5,  1 h,  etc.  are  those  of  the  expansion  of 

|(1  ♦  It  *)*. 

The  subscripts  p  and  e,  suggesting  pole  and  equator,  do  not 
mean  that  we  are  actually  computing  for  the  pole  or  equator,  but  merely 
that  the  notation  is  suggested  by  the  similarity  in  form  to  sinaa  in 
(81)  and  tan2 a  in  (80), 

Series  Expansions  of  and  qi 


q2(a)  -  |  [(3  cot2a  ♦  l)  a  -  3  cota) 


we  find 


■Jfai  A  .  q^  m  3  cot2a  +  2  -  3  a  cot  a  cosec2a. 


These  expressions  are  very  troublesome  to  compute  when  a  is  small. 
He  use  by  preference  the  expansions  in  series  with 


tan  a  ■  c 


s(a)  -  «*  -  A-  «*  +  IsL  -  •  •  •  1 - gS£  +  •  .  • 

8  3-5  5.7  7.9  (2n+l  )(2n+  5)  * 


ql(u)  -  6(^5  -  j-j  *  ~  -  *  •  •  * 


If  we  take  out  the  first  tarn  as  a  factor,  we  have 


qa(«)  •  |jc*  -  f  **  ♦  j  <4  "  ff  •* 


♦  •  •  • 


♦  •  •  • 
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The  quantities  in  square  brackets  may  be  expressed  as  hypergeotnetric 


series  F(o,  Pf  Y;  x)  where 


FK  ti  rs  x)  -  1  ♦  ^  - ** 

'  ’  *’  T’  J  i  •  r  i  *  2  •  r  *  (r+ 1) 


a+2)  • 


1  •  2  .  3  •  Y  •  (Y+l)  (Y+2) 


«,(«)  -  ^  F(f,  2,  |j  -  «’) 


lj(«)  "  |«*F(1,  |»  \i  -  «’) 


Evidently  a  and  p  may  be  interchanged  at  pleasure. 

The  introduction  of  the  hypergeometric  function  enables  us  to 
derive  for  qj  and  qi  other  series  that  may  be  useful  in  checking 
numerical  computation.  The  new  series  involve  quantities  familiar  to 
geodesists.  In  terms  of  the  semi-axes  a  and  b  of  the  meridian  ellipse, 
t  is  evidently  the  second  eccentricity,  or 

c*  -  a  y-b-  -  tan*  a.  (90) 

b 

The  ordinary  eccentricity  e  is  given  by 

e*  ■  -  u"b  -  sin*o  .  (91) 

a 

The  flattening  (ellipticity)  f  is  given  by 

f  »  a- -  -  2  sin*  %  .  (92) 

a  £ 
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There  is  still  another  quantity  n,  which  seems  to  have  no 


special  name;  it  is  given  by 


a  -  b  .  2  a 

n  *  -  *  tan*  -  . 

a  +  b  2 


A  general  transformation  (A)  given  in  elementary  treatments 
of  the  hypergeome trie  function  is 

F  K  P}  Y}  x)  =  (l-x)-a  F  (a,  Y-t}  Yf  (A) 

x  -  1 

Two  rather  special  transformations  are  given  by  Kunmer  11 / . 

F  K  a  ♦  Y;  x)  -  (--  +-^^-f2a  F  (2a,  2a-Y-l;  Y,  <B) 


F  (a,  a  +  i;  Y;  x)  -  (1  -  x)“a  F  (2a,  2Y-2a-l;  Y;  (c; 

2vl-x 

Transforms  (A),  (B),  and  (C)  give  us  respectively  power  series 
2  1 

in  e  ,  n,  and  — f  with  numerical  coefficients  as  follows : 


16  fr  1  2  a  50  25  4 

~  n  1  -  n  +  —  n  -  n3  +  — —  n  - 

15  L  7  7  231  11*3 

2_  3  T,  +  9  8  f 2  +  200  f3  600  . 

15  6  L  7  f  +  7  231  1001  ^  * 


1X/  E.  E.  KuBaer.  "liber  die  hypergeometrische  Reihe 

1  ♦  ZlL  x  ♦  Mill  +  "Ca»D  (g*2)  P( 


I +1MP  +21  * 


1*T  1  •  2  •  Y •  (Y+l)  1.  2.  3  Y  (r+D  (r+2) 


♦  •  «  •  w 


(Crelle's)  Journal  fur  die  reine  und  angewandte  MatheaatUu  1036.  Vol  15, 
39-8 3.  Hummer  gives  a  great  number  of  transformations.  These  used  here 
are  numbers  1|2  and  1*7  on  page  77. 
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*  |L  e»  ♦  M.  e»  ♦  .  .  . 

231  3003  J  ’ 


o'  «Se*r  l  +  -e2+-^-a4  +  e6  ♦ 

r  L1  ?  21  231  3003 

.  4  n  P  1  ♦  —  n  -  “  +  A.  n3 - n4 

5  L  7  21  231  3003 


2  ,  r  8  20  . 

5  e  L1  +  7f  +  af 


2  +  l60  »a  +  200  ». 

231  U29 


*  .  .  . 


The  series  in  n  are  the  most  advantageous  for  computation, 
since  the  signs  alternate  and  for  small  values  of  a  we  have  roughly 
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Table  I 


1 

f 

ao 

qa(“o> 

9a(a0) 

296.0 

4°  42 •  39.6* 

.0000741879 

.00270858 

.1 

36.8 

1502 

767 

.2 

33.9 

1126 

675 

.3 

31.0 

0751 

584 

.4 

28.2 

0375 

492 

296.5 

25.3 

0000 

401 

.6 

22.5 

.0000739626 

309 

.7 

19.6 

9251 

218 

.8 

16.8 

8877 

127 

•9 

13.9 

8504 

036 

297.0 

11.0 

8130 

.00269944 

.1 

8.2 

7757 

853 

.2 

5.3 

7384 

762 

•  3 

2.5 

7012 

671 

.4 

4°  41*  59.?" 

6640 

581 

297.5 

56.8 

6268 

490 

.6 

54.0 

5897 

399 

.7 

51.1 

5525 

308 

.8 

48.3 

5154 

218 

.9 

45.4 

4784 

127 

298.0 

42.6 

4414 

037 

.1 

39.8 

4044 

.00268946 

.2 

36.9 

3674 

856 

•  3 

34.1 

3305 

765 

.4 

31.3 

2936 

675 

298.5 

28.4 

2567 

585 

Table  la 

Spheroid 

ao 

9a(°ho) 

q'a(ao) 

Clarke,  1880 

4° 

43'  52.8" 

.0000751524 

.00273203 

Clarke,  1866 

4° 

43'  09.0" 

745739 

271798 

Bessel,  1841 

4° 

41*  10.0" 

730167 

267998 
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Table  II 


Values  of  a 


Ar<  0°  45°  90° 


0 

4° 

42' 

11.0" 

4° 

42' 

11.0" 

4° 

42' 

11.0" 

10 

4° 

41' 

44.5" 

4° 

41' 

44.5" 

4° 

41' 

44.6" 

20 

4° 

41' 

18. 0" 

4° 

41' 

18.1" 

4° 

41 ' 

18.2" 

30 

4° 

40' 

51.6" 

4° 

40 ' 

51*7" 

4° 

40' 

51.9" 

40 

4° 

40' 

25.3" 

4° 

40' 

25.5" 

4° 

40' 

25.6" 

50 

4° 

39' 

59.1" 

4® 

39’ 

59.3" 

4° 

39' 

59-5" 

60 

4° 

39' 

32.9" 

4° 

39' 

33.2" 

46 

39' 

33.4" 

70 

4° 

39' 

06.8" 

4° 

39' 

07.2" 

4° 

39' 

07.5" 

80 

4C 

38' 

40.9" 

4° 

38' 

41.2" 

4® 

38' 

41.6" 

90 

4° 

38' 

15. 0" 

4° 

38' 

15.3" 

4° 

38' 

15.7" 

100 

4° 

37' 

49. 1" 

4° 

37’ 

49 . 6" 

4° 

37' 

50. 0" 

Values 

of  q2(a) 

0° 

45° 

90° 

0 

.0000738130 

.0000738130 

.0000738130 

10 

4659 

4665 

4671 

20 

1210 

1222 

1233 

30 

727783 

727800 

727817 

40 

4377 

4400 

4422 

50 

0992 

1020 

1048 

60 

717628 

717662 

717695 

70 

4285 

4324 

4363 

80 

0963 

1007 

1051 

90 

707661 

707711 

707760 

100 

43 80 

4435 

4489 

Values 

of  qg(a) 

0° 

45° 

90* 

0 

.00269944 

.00269944 

. 00269944 

10 

9097 

9098 

9100 

20 

8253 

8256 

8259 

30 

7413 

7418 

7422 

40 

6578 

6583 

6589 

50 

5746 

5753 

5760 

60 

4918 

4926 

4934 

70 

4094 

4103 

4113 

80 

3273 

3284 

3295 

90 

2457 

2469 

2481 

100 

1644 

I658 

1671 

The  values  of  fir  are  given 

in  kilometers 

.  Note  that  at 

equator  and  pole  Ar  represents  an  altitude,  but  at  geo- 

centric 

latitude  45®  it  differs  from  an 

altitude,  which 

incidentally  might  be  variously  defined, 

by  a  quantity 

of  the 

order  of  If3&r> 
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